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Abstract: It is old folklore that the violation of Leibniz rule on a lattice is an obstruction 
for constructing a lattice supersymmetric model. While it is still true for full supersymme- 
try, we show that a slightly modified form of the Leibniz rule, which we call cyclic Leibniz 
rule (CLR), is actually a criterion for the existence of partial lattice supersymmetry. In 
one dimension, we find sets of lattice difference operator and field multiplication smeared 
over lattice which satisfy the CLR under some natural assumptions such as translational 
invariance and locality. Thereby we construct a model of supersymmetric lattice quantum 
mechanics without spoiling locality. The CLR relation is coincident with the condition 
that the vanishing of the so-called surface term in the construction by lattice Nicolai map. 
We can construct superfield formalism with arbitrary superpotential. This also enables 
us to apply safely a localization technique to our model, because the kinetic term and 
the interaction terms of our model are independently invariant under the supersymmetry 
transformation. A preliminary attempt in finding a solution for the higher dimensional 
case is also discussed. 
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1. Introduction 

Supersymmetry (SUSY) is not only a candidate for physics beyond the standard model 
but also a possible framework (or a part of it) with definite ultraviolet behaviors in 
quantum field theories. If we can construct supersymmetric theories nonperturbatively, 
we shall obtain deeper understandings of vacuum property and other nonperturbative 
behaviors [[[], [2|, |3|. Although the SUSY algebra is generally incompatible with a discrete 
lattice space, there is no obstruction to formulate a free SUSY theory on lattice, if we 
replace an infinitesimal translation with a lattice translation. 

Our real difficulty to overcome is Leibniz rule (LR) on lattice in interacting SUSY 
theories. Unfortunately, the no-go theorem for the rule on lattice has been proved in 
our previous papers j|, [|. The essential point of the proof is the incompatibility of the 
translational invariance, the locality and the LR on lattice. Under a certain circumstances, 
SLAC type difference operator satisfies the LR on lattice^, ||| but it is non-local. 

As another approach, we could utilize a Nicolai map for constructing a lattice SUSY 
theory^]. For supersymmetric quantum mechanics, a lattice model with a Nicolai map 
was proposed in Ref.[||], 

SN = E \ (J A< ^ " W M) 2 + E it™ Umn + V>n , (1.1) 
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where A mn is a difference operator and W n {(j)) is an arbitrary function of 4> m . The lattice 
action possesses an exact sup er symmetry fi |l0| but contains the unconventional term 

-^(A0) n W„(0). (1.2) 
n 

It should be noticed that the above term is absent from the action of supersymmetric 
quantum mechanics and would vanish in the continuum as a surface term. It cannot be, 
however, written into the form of total divergences due to the lack of the Leibniz rule on 
lattice. The author [ 11 1 has discussed, in a topological point of view, an off-shell action 

Sb = E l^ 2 + E ZjB «(( A ^« - w m) + E^™ f A ™ + ^r^-) ^ > (°) 



which reduces to (LI) after eliminating the auxiliary field B n . The above lattice action 



was rederived in a superfield formulation [ 12 1. Although the lattice action has an exact 



fermionic symmetry, it may be regarded as a lattice version of the stochastic action [13 1 
rather than the standard supersymmetric quantum mechanics, 

*»» - / dt [ ) 2 + \ f " + <™™ + * a + w ™ 

with a Euclidean time t. Hence, it will be worth looking for a lattice formulation which 
leads to a lattice version of Ssqm and also does not rely on the existence of local Nicolai 
maps without unwanted terms like ( |1.2| ). 

In this paper, we propose a novel criterion for the lattice SUSY realization which we 
call cyclic Leibniz rule (CLR), that is not a genuine LR on lattice but still a natural lattice 
analog of the Leibniz rule in the continuum theories in the sense that it is reduced to the 
corresponding relation in the continuum limit. We construct sets of the lattice difference 
operator and the associated symmetric field product which satisfy the CLR. With these set 
of difference operator and field product, we formulate the concrete example of lattice SUSY 
quantum mechanics. Although only a part of SUSY could be realized exactly, yet they have 
several advantages over the existing approaches. One of the advantages is that our lattice 
model is free from the "surface" term ( |1.2j ) and is just a lattice analog of the standard 
supersymmetric action (|l-4)- This property with CLR allows us to construct a superfield 



formulation, whose realization is different from that given in Ref.[12|. In the action of 
our model the kinetic term and the interaction terms are separately invariant under SUSY 
transformation. Because of this property we can apply a localization technique to obtain 
some exact results in our model. Since our lattice formulation mimics the continuum 
supersymmetric quantum mechanics and does not rely on the existence of local Nicolai 
maps, it may shed a new light on the construction of higher-dimensional supersymmetric 
lattice models. 

In section 2, we propose the cyclic Leibniz rule (CLR) as a criterion for the SUSY by 
using a concrete example, for definiteness, of D = 1 J\f = 2 Wess-Zumino model on lattice. 
The lattice difference opertors and associated symmetric product is found as a solution 
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to the CLR. Based on these, in section 3, we construct exactly (half-)supersymmetric 
quantum mechanics on lattice with arbitrary interactions. The Witten index is exactly 
calculated by the Nicolai map and the localization technique. Multicomponent cases are 
discussed in section 4 and we find a local solution which is a natural generalization of 
one-component case. In section 5, we give an attempt to a higher-dimensional case letting 
spinor index aside, while we so far could not find any suitable solution of the CLR for 
euclidean field theories. Section 6 is devoted to conclusion and discussions. In appendies, 
we find the general solution of CLR with locality, the properties of the associated product, 
and explicit multi-field product examples based on CLR. 

2. A criterion for lattice supersymmetry 

In this section, we begin by recapitulating the difficulty in constructing interacting super- 
symmetric theories on lattice where the important key concept is the Leibniz rule. Our 
proposal to avoid the difficulty will appear subsequently in this section. 

For concreteness, we consider D = 1 model, i.e. supersymmetric quantum mechanics, 
which consists of a multiplet 

(0„,Vn,^n) (2.1) 

where is a real scalar, F is a real auxiliary field and ip is a complex fermion. Index n 
stands for a lattice site, na is a coordinate but hereafter we set lattice constant a = 1 for 
simplicity. The M = 2 SUSY transformation is described as 

S(j) n = etp n - etp n , 
S^ n = e(i(A<f>) n + F n ), 
5tp n = e(-i(A(f)) n + F n ), 

5F n = -ei(A$) n - ez(AV>)n, (2.2) 
where a difference operator A is defined as 

(Av4) n = ^ ] A wm ^4 m , ^ ] A nm = 0, (2-3) 

m m 

for arbitray lattice field A n . The second relation states that the difference of a constant 
field vanishes. The algebra follows from the above definition: 

[Si,62}(4>n,^n,F n ) =2i(€ie 2 -e 2 ei)((A^) n ,(A^)„,(AF) n ). (2.4) 
The kinetic terms Sq and the mass term S m of the free action 

S = (\(^4>)l + *(A^)„ + l -F 2 n ), (2.5) 

n 

S m = im^2 ( F n4>n + ^n^n) (2.6) 
n 



-3- 



are invariant under ( p.2|) provided 

^nm — A mn (2'7) 

which is a generalized symmetric difference operator. 

The obstruction for SUSY realization on lattice appears in the interaction parts. For 
an illustration we consider here cubic interaction such as 

s int = e| (M fmn F e 4> m 4> n + 2N £mn (f) e ip m ifj n ), (2.8) 

where g is a coupling constant, and Mg mn and Ni mn are coefficients for fields products. 
We can impose Mi mn = Mi nm without the loss of generality. Under fl2.2|) , the interaction 
terms are transformed as follows, 

Imn 

-2e(M£ mn - N m £ n )F e ip m (f) n + 2e(Me mn - N nm i)F e ^ m (j) n 

+ 2A^ mn ^^ m (A0) n ) I . (2-9) 



In order to accomplish the invariance under (|2.2| ), we must require 



^2(A ke M kmn + N nik A km + N mek A kn ) = 0. (2.10) 
k 

for e invariance and 

Me mn = N m £ n , 

y^j(A k eM kmn + N nk iA km + N mk iA kn ) = 0, (2-11) 

k 

for e invariance. These conditions can be combined into a simpler form 

Me mn = M£ nrn = M nrn £ = N£ mn = Ninm = N nm £, (2-12) 

y^j(A k eM kmn + M£ kn A km + M^ m fcAfc„) = 0. (2-13) 

k 

With (|2.7|), the relation ( |2.13| ) implies a Leibniz rule on lattice 

Aik X] M kmn^m^n = ^ M ekn (A(j)) k 1p n + ^ M £mk (j> m (Alp) k , (2.14) 

where ^ are any lattice fields. 
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Before going into the further analysis, we describe the translational invariant and local 
difference operator and product-coefficients as holomorphic functions jl. From transla- 
tional invariance, these objects depend only on the difference of site indices, so we denote 
conveniently 

A(k) = A rnm _ k , M(kJ) = M mm _ km _ e . (2.15) 

By transforming them to the momentum representation (w— expression), the complex func- 
tions 

oo oo 

A(w) = w n A(n), M(w,z)= ^ w m z n M(m,n) = M{z,w) (2.16) 

n=— oo m,n=— oo 



are defined. Substituting w = exp(ipa), z = exp(iga), the expression ( p. 16 ) exactly corre- 



sponds to its Fourier expansion and p, q become momentum variables when p, q are real. 
The locality of the difference operator is exactly equivalent to the holomorphism of the 
^-expression, and the holomorphic domain T> is an annulus in the complex w-plane which 
includes a unit circle around the origin. For the product M(w, z), the domain is T> <g> T> as 
a two- variables function. Namely, we proved that the locality, i.e. long distance behaviors 



I A m n 



< Cexp(-K \m-n\), \M tmn \ < C x exp {-K x \l - m\ - K 2 \£ - n\) (2.17) 

with C, Ci, Kq, Ki, K2 > are equivalent to the holomorphism of A(w) and M(w, 
From the second relation of fl2.3|), the correspondence A has the property, 

A(l) = 0. (2.18) 

Turning the relation ( |2.13| ) into the ^-expression, we can find 

M(w, z)A(-) + M(w, z)A(w) + M(w, z)A(z) = 0. (2.19) 
wz 



The only holomorphic solution around w = z = 1 of ( 2.19j ) is A(w) ~ logw i.e. a SLAC 



type difference operator besides a trivial solution with M(w, z) = 0. The SLAC-type 
operator is, however, not holomorphic in the annulus T>. Thus there is no set of A and 
M(^ 0) with translational invariance and locality that satisfies ( |2.13| ). This is an essence 
of the no-go theorem for the Leibniz rule on lattice proved in our previous work|3j]. So the 
lack of lattice Leibniz rule is an obstruction for lattice supersymmetry with interaction. 

To overcome the above problem, we demand only a subset of supersymmetry. For 
example, if we impose the invariance only for e in (^]^), then we have the half of the 
condition, i.e. (|2T0|) without ( [Till) . Thus, instead of (|1|) and ( ^T3|) , we obtain 



Me mn — Mi nra — N nm £ — N mn £, (2.20) 
y^(AfclMfc mn + MknlAkm + MkirnAkn) = 0. (2-21) 



We call ( |2.21| ) as cyclic Leibniz rule (CLR). This relation cannot be factorized, like ( 2.19| ), 



into a separate condition for a difference operator, but must be solved by the combination 
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of a difference operator and a product. It is worth noting that the symmetry in the product 
Mkim of ( 2.12j ) is different from that of ( 2.20| ) ; the former is totally-symmetric for the three 
indices and the latter is only symmetric for the second and the third indices. Although the 
conventional LR ( 2.1 3| ) seems equivalent to the CLR ( 2.21| ) in the naive continuum limit 
because the limit implies M k £ m —¥ 5 k n5 km , they clearly have different property for the finite 



lattice constant. In contrast to (2.14), the CLR can be expressed as 



X] A ^ X] M kmn<Prn^n = ^ M kn£ (A<p) k 1p n + ^ M klm (\) m ( Alp) k , (2.22) 



n,k 



k,m 



whose indices are indeed summed in a different way from ( |2.14| ). Although the CLR is 
found in constructing a specific model with cubic interaction, it will turn out to be very 
useful for the analysis of more general super symmetric theories in later sections. Instead 
of the usual Leibniz rule fl2.13|) , we can adopt the CLR as the key rule for the exact lattice 
supersymmetry. As seen later, we can construct a superfield formalism and find a Nicolai 
map with vanishing surface terms in interacting cases. This framework with the exact 
symmetry are a real advantage of our approach . 

In w-represenation, the cyclic Leibniz rule ( |2.21| ) can be expressed as 



M(w, z)A(—) + M(z, —)A(w) + M(—,w)A(z) = 0. 

wz wz wz 



(2.23) 



We give here an example solution for (2.23) keeping translational invariance and locality 



Amn. 



(25, 



Lm-l<- 



',ri-l + ^,n-l^,m+l + 



(5m,n— 1 <^m,ra+l)- 



M(w, z) = - (2wz + wz 1 + zw 1 + 2{wz) 1 ), A(ir) 

This corresponds to a real space expression 
1 



w — w 



',n+l°e,m-l 



25, 



e,m+ioe, n +i), 



(2.24) 



(2.25) 



This simple example has doublers which can be resolved by a super symmetric Wilson 
term or a super symmetric overlap-type mass. We shall see these terms explicitly in the 
subsection 3.1. More generic other solutions to (|2.23|) are given in Appendix [Al 



3. Analysis of supersymmetric lattice quantum mechanics 

In the previous section, we have found the important relation, the CLR. Its solution con- 
sists of a local difference operator and a symmetric field product smeared over the lattice 
sites keeping locality. In this section, we explicitly construct general supersymmetric quan- 
tum mechanics using the CLR as a guiding principle. And we will analyze the model by 
localization technique. 
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3.1 Construction of supersymmetric action 

We saw in the previous section that a cubic action can be half-super symmetric if the prod- 
uct satisfies the CLR. This can be immediately extended to other order terms including 
quadratic, quartic or higher. An important point is that the CLR for each order guarantees 
the invariance of each order term independently. In other words, the coupling constant for 
each order term can be taken independently. This naturally leads us to the superfield 
formalism with half-supersymmetry to construct generic action. Also we are able to uti- 
lize this property to make an exact calculation with localization technique which will be 
discussed in the next subsection. 

Let us first define a 2-body product of lattice fields using the coefficient Mi mn 

{cf>, if,} t = ^ M £mn <p m ^ n = (-1) ^{V, (3-1) 

mn 

where and are the Grassmann-parity for <f> and ip. We refer this product as a smeared 
symmetric product, because this is a symmetric product smeared over lattice sites. The 
locality of the product is guaranteed by the holomorphic property of its ^-expression 
M(z,w). We also define an inner product of two lattice fields 

(<£,V)=l>nV>n. (3.2) 

n 

Then a simple 3-body summed product of lattice fields is expressed as 

(3-3) 

Imn 

With these notations, as well as (|2.7j), fl2.2l|) , the CLR can be rewritten into a compact 
form 

{V, x}) + (AV, {X, </>}) + (Ax, {</>, V>}) = 0. (3.4) 
In the case of three same fields, ([^|) is reduced to 

(A<M<M}) = 0. (3-5) 

The expression Q3.5| ) is equivalent to ([O]) if we take the symmetry of the product M into 
account. These simple expressions are useful in further SUSY analysis. 

Now we construct various terms in the action using superfield formalism. As stated 
before, the kinetic term is fully supersymmetric with e and e, while the interaction term is 
half-supersymmetric with, say e. Therefore we can use N = 2 superspace for the former, 
but must use M = 1 for the latter. 

An M =2 supermultiplet for a superspace (n, 9, 6) 

3 n (0,0) = * n (0)-0tf„(0) (3.6) 
is decomposed into two M =1 multiplets, 

= ^n + (3.7) 
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which is Grassmann even and 

* n (0)=il>n + eF n (3.8) 

which is Grassmann odd. 

A kinetic term of the action can be written as 

So = \J Sd9(DE(9, 9), DE(9, 5)), (3.9) 

where supercovariant difference operators are defined as 

. d - . d 

Dfnn = ^TT^^mn 9A mn , D mn = 2— =(5 mrl 9A mn . (3.10) 
69 69 

To formulate mass term and Wilson term supersymmetrically, we consider a generic 
bilinear form 

(<f>,Ghl>) ^YtfaGwPt, (3.11) 

kl 

where G mn is a certain operator on a lattice. The translational invariance implies that G 
is a function of the difference of its indices G mn = G{m — n), and satisfies AG = GA with 
translationally-invariant difference operator A. Thus we have 

(A0,Gty) = -{</>, GAj>), (3.12) 

with A T = - A. Then 

S 2 = i J d9(V(9),G<S>(9)) (3.13) 

gives a super symmetric quadratic action invariant under e-transformation. For example, a 
simple mass term with Wilson term is given by choosing 

G(w) = Y J W n G{n) = im + ir 2 ~ W ~ W 1 (3.14) 
n ~ 

where m is a bare mass and r is a standard Wilson term parameter. For more generic 
doubler-suppressing term, we can write 

G(w) = im + iH(w,M), (3.15) 

where H(w, M) is a holomorphic function in T> with a parameter M > and corre- 
sponds to w-expression of a translationally invariant operator H mn (M), i.e. H(w,M) = 
J2m wm Hee+m{M). We further impose that H is symmetric, i.e. H T = H or H(w) = 
H(l/w), throughout this paper. If we need an overlap-type action, we choose H to have 
the poles at w = 0, oo and the branch points at w = —1/M, —M with M > 1, as well as the 
property H(w = 1,M) = 0. The explicit form for the overlap-type term||i"4|, |l5| 1 is given as 

H(w,M) = 1 (M + H+Hli) = i _ (M+jgifZi) 

M+l I ( w-w-i- yi | | ^+t»-^ 2 v / (M + w)(M + w" 1 )' 
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Note that, together with an explicit solution for the CLR 

*, . (M + ^tw-w- 1 ) 

A(w) = , , (3.17) 

2y/{M + w){M + W- 1 ) V ' 

M(w, z) = — — - ry/(M + w)(M + W - l )(2wz + wz' 1 + zw~ x + 2(wzY x \ , (3.18) 

6{M + 1) V / 

the above H satisfies an analog of the Ginsparg- Wilson relation) 16, 17], 

A(w)A( W - 1 ) + H{w, M)H(w~ 1 1 M) = (M + 1) (h(w, M) + H^" 1 , M)) , (3.19) 



where the Dirac operator is given by DDirac(w) = iA(w) + iH(w,M). 

The quadratic terms including mass and doubler-suppressing term and the cubic in- 
teraction term of the supersymmetric action are now obtained as 

s 2 + s int = i J de((*{o),GQ(e)) + g2(*(e),wo)Mo)}j) ( 3 - 20 ) 

with a coupling constant g 2 (here subscript 2 stands for 2-body product). The action ( 3.2C| ) 
is clearly invariant under SUSY transformation 

8cj) n = e ^ n , 

5i/j n = e(i(A<f>) n + F n ), 
5ip n = 0, 

8F n = -ie{Ai)) n . (3.21) 
This transformation can be realized on the superfield as follows, 

6E n (9,9) =J2 e {9Q S nm + i0A nm jE m (e,e) = e(— $ n (0) + 9— ¥ n (0) + i9(A$) n (9)), 

m 

(3.22) 

or equivalently on the M = 1 superfields 

<5* n (0) = ^n{9) + ie(A$) n (9). (3.23) 

Thus the total action 



St — So + S 2 + Sj 



mt 



= ~(A&A0) + i$,A^) + i(F,F) 

+ (F, (fm + iH)<j>) + (*m + iH» + 132 (F, {0, 0}) - 2i 52 (V, {0, ^}) (3.24) 
is also invariant under ( |3.22j ) or ( |3.23j ), 

SS t = 0, (3.25) 
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due to the relation ( 3.1 2| ) and the CLR ( |3.5[) . It is worth emphasizing here that our for- 



mulation keeps locality; the holomorphic property of the ^-expression of various operators 
and product guarantees the locality. 

To extend to multi-body supersymmetric interactions, we define the smeared iV-body 
products for a bosonic superfield 

;mi,...,mjv 

V'"^, (3.26) 

mi,...,mjv 

where mi, .., m^v are totally symmetric indices (we have been omitting semicolon for N = 2, 

i.6. M-n\m\,m,2 — ^ntnim2 ■) 

If this product satisfies the iV-body CLR 

(A*, {*,... ,*})=0, (3.27) 

then the action 

J dOig N (*,{*,...,*}) (3.28) 



is invariant under the supersymmetry transformations ( |3.23| ) in exactly same way as the 
2-body case. It should be noticed that the i?-symmetry 

e -> e ia e , 
$ -»• $, 

* -»• e ia * , (3.29) 
allows only the type of interaction terms ( 3.2Sj ), which consists of <!> and VP and is necessary 



for constructing supersymmetric 2iV-body interactions for scalar fields with a coupling 
constant g^. Also note that iV-body product satisfying the CLR can be constructed in 
terms of 2-body products which will be discussed in Appendix B. 

Before closing this subsection, it may be helpful to understand the exact SUSY invari- 
ance in an alternative way. We write the total action St as a sum of kinetic term So, bare 
mass term S m , doubler-suppressing term S^ s and interaction term S{ n t 

St = 5*0 + S m + Sds + Shit- (3.30) 
Then each term can be expressed as a supersymmetric transform of something 

So = V^HA^l + F)), 
S m = im5'(ip,(j)), 
S ds = id'^H^), 

S int = S 3 + S 4 + --- = ig 2 6'(ip, {</>, 0}) + ig 3 6'(ip, {0, 0, 0}) + • • • , (3.31) 

with coupling constants <72> <73> • " "• Here we denote the transformation without parameter 
by 5' , i.e. 5 = ed' , or explicitly 



'71 • 



8'4> n = 1p n 

<5>„ = »(A0) n + F n , 

5'i>n = 0, 

6'F n = -i(A^) n . (3.32) 
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Due to the nilpotencey of the transformation 5' 2 = 0, the invariance of each term is 
manifest, 

5' So = 5'S m = 8'S ds = 5'S int = 0. (3.33) 
3.2 Localization technique and calculation of Witten index 

With the action ( fOoD , we can find a Nicolai map for our system, and apply the localiza- 
tion technique fl~8|, |l9[l to calculate the Witten index (the partition function with periodic 
boundary condition for all variables). 

Before going into the detail, it would be useful to compare with the previous attempts 
in the existing literature. In the references |^, |l2| the partition function has been calculated 
in the ultra local limit, where the iV-site partition function reduces to N copies of the 
one-site one, i.e. Zn = (Z\) . This result indicates that every degrees of freedom may 
contribute to the Witten index that is quite different from the usual understanding of its 
topological nature. Also it depends on whether the number of lattice sites is even or odd 
for Z\ = — 1 case of which we could not understand the validity. On the other hand, our 
lattice formulation is suitable to use the localization technique as will be seen below. It 
will turn out that the Witten index is determined by only zero modes and all fluctuation 
modes are irrelevant in it. 

Now let us begin by defining the Nicolai map 

in = ~(A(p) n + m(j) n + {H<j)) n + g 2 {(j>, <P}n + 93{<f>, <P, <P}n H • (3.34) 

With this, the Witten index can be calculated as 

Z = j D<j)D^D^DFe~ St = j D£D^D^\ detr/]" 1 e~ s ™ §^«-*Em VroW'fc 

= J D£\detri\- 1 detr] e _s »^" $ ", (3.35) 
where all variables follow a periodic boundary condition and 

= ^— = + m5 ik + (H) ek + g 2 -^—{4>, (f\k + 93Trr{4>, <t>, 4>}k H ■ (3.36) 

o<p£ ■ ocpe o(p£ 

The prefactor of ( |3.35| ) 

| detT/l -1 det?7, (3.37) 

is not generally a constant but det r/ can vanish for certain cases. Also the path of the 
integral ( 3.35| ) is not definite. Therefore we must carry out more detailed analysis. 



In order to control the infrared behavior, we restrict a size of the system to finite N. 



Also we assume the highest power of the potential in (3.34) is p, for definiteness. We rescale 



our kinetic and doubler-suppressing term with parameters t and ids respectively, 

So — > tSo, Sds — > tdsSds- (3.38) 
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Note that t ds looks like a Wilson term parameter for Wilson fermion case. Then we integrate 
out F n and obtain a i-scaled partition function, 

Z(t,t ds ) = J Dcf)D^D^DF e - s{t ' tds) 

= [ Dcf)D4'D4'DF e~ tSo ~ tdsSds ~ Sm ~ Sint 
t -N/2 i D( p D ^ D ^ 



x e 



-(^+*L^ 1 ^ s +t- 1 ^ +int +t ds t- 1 S c 6 ross )-(iS i (+t ds 5£+5^ +int ) 



(3.39) 



Here we denote 5q, 5^,5j s and Sf nt as bosonic (without fermion) parts of kinetic, mass, 
doubler-suppressing and interaction terms, respectively. S[ , Sm, S[ s and S[ nt are fermionic 
(with fermion) counterparts. We summarize each part of our action in Table |]. In addition, 

Table 1: Bosonic and fermionic parts of the action after integrating of F n 



ob 


ob 


ob 
m+int 


±(A0,A0) 




\(m<j) + g 2 {4>, 4>] H ,m4> + g 2 {4>, <ft} H ) 













im(Tp, ip) - i(ip, 2g 2 {(p, ^} H ) 



the cross terms S^ross for doubler-suppressing term and mass plus interaction terms are 
defined as 

5 c fe ross = {H<f> , m<P + g 2 {<P, 0} + • • •)• (3.40) 



It is clear that this Z(t, t ds ) is t-independent, because before -F n -integration one can 



show 



§- t z(t,t ds 



J DcpDipDipDF S e-^o-Sm-tdsSds-Smt 

D(j)D^D^DF5'(Xe~ tSo ~ Sm ' tdsSds ~ Sint ) = 0, 



(3.41) 







at. 



ds 



■Z(t, t ds 



J DcpDipDipDF S ds e - ts o-s m -t ds s ds -s int 

D(j)D$D7pDF5'(Ye~ tSo - Sm - tdsSds - Sint ) = 0, (3.42) 



where X = ^ yp, —iA<J) + FJ , Y = yp, Hcpj . Here we used the invariance of the path- 
integral measure under 5' and 5' (S(t,t ds )) = 0. 
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In the following calculations of the Witten index, we take two typical limiting cases. 
One is analogous to the continuum theory in which important modes are only zero-modes. 
The second case is peculiar to lattice theories, which has doubling modes in addition to 
zero- modes. 1 Although the two cases are physically and mathematically different, after 
using localization technique, we can obtain the exactly same results. We first decompose 
each field into zero mode (script 0), doubler mode (script d) and the other modes: 

6n (-1)" , 

H 7=-<Pd + a n , 



N VN 

ipn = -j= H TrT^d + Xn, 

VN VN 



7 ^0 . ("I)" T I — /« ,q\ 

V>n = -7= H nFri>d + Xn (3.43) 

V N V N 

Then we rescale each part with t or t^ s in such a way that the only important mode survives 
in the infinitely scaling limit. 

Case (i) We make the scaling of each field as 

4>d + ~^ a rn 

• •» - — t= + V — ipd + ~qXn, 

*" = # + v^*< + 7P' (3 ' 44) 

then we take t,t^s —> oo with a fixed t/t^g- This case is similar to the continuum theory, 
because the real zero mode (not doubler) is only t, ids-unscaled. A scalar zero mode <j>Q 
contributes to t~ 1 S^ n+int + £^ +int and a scalar doubler mode <f>d contributes to 



<f>0 


+ - 


(-i: 


P 


Vn 




ids 


^0 


+ - 


(-i; 


P 


Vn 


\/Ny 


ids 


^0 


+ - 


(-i; 


p 


VN 




ids 



= %E(E^-(- 1 ) n ^) 2 = ^(-l)) 2 ^, (3-45) 



t as tN ^ \^ ' T V t 

m n 



where H(w) = ^2 n w n H(n) for translational invariant H mn = H{m — n). On the other 
haj 
to 

t.^sl = % 

N 



hand, the fermion zero modes contribute to S m+[nt , and a fermion doubler mode contributes 
ids^L = 4 ^2(-l) m ^dH m n(-l) n ^d = iH(-l)M d . (3.46) 



Any non-zero modes can be integrated out in the limit t, ids ~~ * 00 with fixed t/t^s, and 
the t factor of the partition function ( |3.39| ) is exactly canceled. (t^ s /t)S^ TOBS terms in the 
action vanish in the limit. 



1 Hereafter we assume the number of sites N is even for simplicity. There is no doubler mode for odd N, 
so that it gives the same result as Case (i) below. 
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Therefore 

Z (t,t ds ) = ^ J d^d<p d d^d^ d di> d e -(^*- 1 ^+*- 1 ^ + »t)-(*ds5 d ^+5^ +int ) 



-^H(-l)Vt f de^rS) e"^ 2 , (3.47) 

1 ) ) « 



where 



t 

* Vt ds (H(-l)) 



^(^o) = m + %M(1, l)0o + ^M(l, 1, + • • • + -y^Mil, • • • , l)^" 1 (3.48) 



with ^-expression of translaionally invariant multi-body product 

M(wi, ■■■ , w n -x) = Yl ■■■Wn-i 1 M t/+mi ,., i+mn _ 1 . (3.49) 



mi,-,mn-i 



The final expression ( |3.47 ) results in Z{t,t ds ) = for p =even and Z(t,t ds ) = ±1 for 



p =odd (its sign depends on the asymptotic behavior of the potential). 
Case (ii) An alternative way of rescaling of fields is given by 

1 jl (-!)".. 1 
- - i> H 7=r4>d + —pa n , 



i _ (-D™ _ i 

V>„ = -7=^0 + ^7=-V>d + -JjXn- (3.50) 
VJV vJV v* 

Then we take the limit t — > oo with a fixed t ds - In this case not only zero mode but also 
doubler mode are t- unsealed. A 0o contributes to 

* ^m+inti ^m+int' (3.51) 

and a 4> d contributes to 

^ds^ •S'dsi t S m+int , td s t S cross , 5'^ l _ ) _ int . (3.52) 
Fermion zero modes tpo^tpo contribute to 

^m+inv (3.53) 
and fermion doubler modes ijjd, tpd contribute to 

^ds'S'ds' ^m+inV (3.54) 
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Any non-zero modes can be integrated out similarly to the case (i) in the limit t 
Therefore Z is evaluated as 
1 



oo. 



Z(t,t d 



x e 



-{t-H 2 , s b , +t- 1 t dB s^,+t- 1 s b ,. t )-(t dB sf +s f , . ,) 

V ds ds 1 as cioss 1 m + int' \ as ds 1 m+int-' 



d^ d \det V (M0, MOT 1 *** v(MO,MO)^ e °'^l (3-55) 



where 



tt(<h,<f>d) = i V^n(-l) n ^ + ^= Vm4„(0 O + (-l)> d ) 

Vt \VN n VN n 

» n i ■ - - n „ 



and 



v(<Po,4>d) 



( d£a(4>a,<t>d) dZojfo^d) ^ 
<9</>o <9^ 

d£d(<j>o,<j>d) dCd{4>o,4>d) 



900 / 
The final expression (EDSBD results in Z(t,t ds ) = for p =even and Z(t, t ds ) 
p =odd. 



(3.56) 



(3.57) 



±1 for 



In summary, these two limiting procedures give exactly the same result, although the 
doubler mode survives in the expression ( |3.55 ) of the limit (ii). The final results reproduce 
known results for the supersymmetric quantum mechanics in the continuum. Also there is 
no subtle dependence on the number of lattice sites. 



4. Multicomponent case 

In order to consider a higher-dimensional extension with spinorial index or a multi-flavor 
extension, we need a framework for multi-component fields. A Wilson term or an overlap- 
type term can be made more realistic in the model with multi-component indices. The 
CLR for the multi-component case, however, is slightly modified and should be treated 
with care. In this section we will explain this in a simple situation where each field has a 
single index which we do not specify with either spinor or flavor at the beginning. 

For the fields i/>£, F®) with single index a, we consider the following supersym- 

metry transformation 

S< = e< 
6^=0 

6F% = -ie(A« 



(4.1) 
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where T is a certain constant matrix on the index with property T = T T ', T 2 = 1. We can 
easily obtain a free action including mass and doubler-suppressing terms invariant under 
the above transformation, 

S = i(A0 a ,A0 a ) + i(^,(rUA^) + i(F a ,F a ) 

= - l -5'(r, (A#T ) + ±5\r, F a ), (4-2) 
S m = im{(j) a ,F a )+im{^ a ,^ a ) = im8'{ct) a ^ a ), (4.3) 
5 ds = i(<j ) a ,H ab F b ) + i(i; a ,H ab ip b ) = i5'(<p a , H ab ip b ), (4.4) 

where H ab is a multi-component extension of H in the previous sections and satisfies 
Yl m { H ab)mn = ^2 m H ab(m - n) = H ab (l) = and ((HYA) ab ) mn = -((HTA) ba ) nm . The 
interaction term 

Sint = i92 ((F a , {0, 4>Y) - 2(r, ^V)) = m S\r, {<t>, 0} a ), (4-5) 
is also invariant under ( |4.1| ) where we used {(ft, ip}f. = M^ n (^Vn and 

^((A#y,{0,0r) = O. (4.6) 

a 

The translational invariance implies that the A and M only depend on the difference 
of site indices, as before, 

A mn = A(m - n), M^ n = M abc (k -m,k- n). (4.7) 

Their w— expressions are 

A[w) = w m A(m), M abc {w, z) = ^ w m z n M abc (m, n). (4.8) 

m mn 

Using the combined quantity 

P abc (w, z) = A(l/ (wz))T ad M dbc (w, z), (4.9) 
the multicomponent CLR is expressed as 

P abc (w,z) + P bca {z,l/{wz)) + P cab (l/(wz),w) = 0, 

(4.10) 

P abc (w, l/w) = 0, P abc (w,z) = P acb (z,w). 

Although solutions for the relation ( |4.10 ) are complicated in general, a simple and special 
solution is given by 

P abc {w,z) = f abc P{w,z), (4.11) 

where f abc is a totally symmetric constant and P(w,z) = A(l/(wz))M(w, z) is a solution 
of single component CLR (2.23). This means 

M abc (w, z) = T ad f dbc M(w, z). (4.12) 
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We comment here that when the above indices a, b, c are somehow interpreted as spino- 
rial indices (like dimensional reduction to 1-dim.), the Wilson and overlap- type [fil| , 15] 
fermions appear in analogous way with higher dimension. The Wilson fermion case is 
realized by a doubler-suppressing term tpHijj where 

H£ n = S ab ir 2 ~ 5m ' n - 1 ~ Sm > n+1 , H(w) ab = r 2 ~ W ~ W ' 1 5 ab . (4.13) 

The Wilson-Dirac operator is written as 

- . . w — w^ 1 ^ 2 — w — u> _1 , . * 

iD WD (w) = r + r , (4.14) 

and the bosonic inverse propagator with the doubler-suppressing term is 



w — w 



-l 



2 /2-w - w- l \2 



Db(w) = {^ J -r^ J = -tsAvdHtsAvdM, (4.15) 

where 75 is a matrix satisfying {75, T} = and 75 = 1. 
For the overlap-type fermion, we take 

AH = W ~™~ 1 K OX) (w), (4.16) 

fabc , 

rabc (n ., ^ _ 2 / 2wz + wz -l + zw -l + 2( wz )- 1 



M abc (w, z) = J - - 2wz + wz~ L + zw~ l + 2{wz)- 1 , (4.17) 

6K OD (w) \ / 

where 

KodH = M + 1 (4.18) 

V ' ^(M + w){M + w~ l ) V 7 

The overlap-Dirac operator is written as 

Won =TA(w) + H(w). (4.19) 



The Ginsparg- Wilson relation [|lq, 17] is written as 

A)D75 + 75 ^0D = DodToDod- (4.20) 

This relation leads us to 

A(w)A(w' 1 ) +H 2 {w) = 2iH(w). (4.21) 
Then, the bosonic inverse propagator with the doubler-suppressing term is expressed as 

DbW = -( (1 -^° dW ) + H\ W ) = 2iH(w). (4.22) 

It is a remarkable fact that the inverse propagator of a bosonic variable is determind by 
the fermion doubler's mass term due to the supersymmetry. 
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5. Higher-dimensional case 

In the previous section we consider a multi-component extension which will give a starting 
point for higher dimensional case. In this section, we consider a contrary case, namely two 
dimensional extension of the CLR without spinor index. Although it is not a real higher 
dimensional extension, we can readily see the difficulty which may exist in the real case. 
Real higher dimensional extension will be discussed elsewhere. 

Let us denote a two dimensional lattice sites by m = (mi, 7712) and complex coordinates 
for two- variables w-expression by w = (w\,W2)- Then we use the following simplified 
notations: 

wz = (wiZi,W2Z 2 ) for w = (wi,W2) and z = (zi,z 2 ), (5.1) 

also we write 

1 = (1,1) and w 1 = (u^.ioj 1 ). (5.2) 
Now we consider the following CLR for any bosonic fields 

(A*, {% 3}) + (A*, {S, $}) + (AH, {<D, vl/}) = 0, (5.3) 

where the two-dimensional difference operator and the symmetric product are defined by 

A = (Ai, A 2 ), {% ~} k ee M kmn*rnZ n . (5.4) 

mn 

Here Ai and A2 operate the first and second site index (mi, 771,2) of the field respectively. 
By imposing the translational invariance, 

A mn = Aim - n), M kmn = M(k - m,k - n), (5.5) 

two holomorphic functions can be defined as two-variables -^-expression 

A(w) = w n A(n), A(l) = 0, (5.6) 
n 

M(w,z) ee ^ w rn z n M(m,n) = M{z,w), (5.7) 
m,n 

where simplified notation w m = w^w™ 2 is used. Then, the w— expression of the CLR 
can be obtained as 

A( W - l z- l )M{w, z) + A(w)M(z, w^z' 1 ) + A(z)M{w~ 1 z- 1 , w) = 0. (5.8) 
Setting z = 1, we obtain 

A(.-) - A(«), (5.9) 

M(l, 10 ) 

and 

A ^ = »>, T^TTT ^ n 0^*' ^^Mw) + M(w, W-')A(z)) . (5.10) 
M(w, z)M(l, (wz)- L ) v / 
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Prom eg, ( |5.10 ), we have 



A(™ 2 ) = 2 y M(w, (wT 2 )Mw). (5.11) 

M{w,w)M(l,(w) l ) 



By the use of the Weierstrass preparation theorem [|20f on multi- variable analytic functions, 
the above difference operator can be parametrized as 

A{w) = K{w) ( Wl - 1 + (u> 2 - l)u(w 2 )) , (5.12) 

where K and u are nonzero and holomorphic near w ~ 1. 

In the following, we shall prove the CLR Q5.8| ) has only a solution A(«;) oc «iiW2 — 1 
with the local difference operator A and the local product coefficient M. Since also in 
this two-dimensional case, the holomorphism is equivalent to the locality, we discuss the 
relation (fT^) based on a theory of multi-variable holomorphic functions. 

Using eqs. (|5l0|), (|5TTl|) repeatedly, we obtain the following recursion formula, 

A(w n ) = nF(w, n)A(w), (5.13) 

where F is a holomorphic function of w and becomes unity in the naive continuum limit, 
w — > 1. From the parametrization ( |5.12| ), the recursion formula is rewritten by 

n i , / n 1W ^ nF(w,n)K(w) t \ 

W 1 -1 + {W 2 ~ 1)U(W 2 ) = K ( w n) \ Wl ~ 1 + ^ ~ l ) U \ W V ) 

n 

= H( Wl -w^(wZ)), (5.14) 

k=l 

where 

w^iw?) =u k n (l- K " 1)«K)) u n = e 2 ^. (5.15) 
From the uniquness of the factorization near the naive continuum limit 1, 

wi - w {n) {w^) = w x - 1 + {w 2 - l)u{w 2 ) (5.16) 



because any considered function is holomorphic around the limit. The equation (5.16) can 
be exactly solved as 

1 - W - u{1) 

u(w 2 ) = 2 —. (5.17) 

w 2 — l 

Imposing the single-valued condition on u(w 2 ), u(l) must be an integer, I. The resulting 
expression is 

A(w) = K(w)(wi - A) = ^p-i^i ~ !)• (5-18) 
w 2 w 2 

Instead of u>i, we can argue the same thing about w 2 of A(w). Since both cases hold, I 
must be unity and we can obtain a general form of two-dimensional difference operator, 

A(w) = K(w)(wi - —) = K ^ W \ w 1 w 2 - 1). (5.19) 

U>2 W 2 
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This general form ( 5. 19f ) is a trivial extension of one-dimensional difference operators and 
the general solution of CLR can be found easily by substitution: w — > W1W2, z — > z\Z2 
in one-dimensional case. This solution is, however, unsuitable for euclidean field theory 
because, with this, the kinetic term is always written as a form of w\W2 or, in the momentum 
representation, p% +P2- 



6. Conclusion and discussion 

When we naively tried to realize the full supersymmetry on lattice, we were required to 
realize a Leibniz rule (LR) on lattice. The no-go theorem, however, prevented us from 
proceeding further. In this sense, the theorem is indeed an origin of the difficulty in lattice 



supersymmetry [^, 21, 22]. On the contrary, realizing exact half-supersymmetry leads us 
a new criterion, i.e. the cyclic Leibniz rule (CLR). The important thing is that the CLR 
actually has solutions which are realized by the set of local difference operator and local 
symmetric field products. Unlike the LR, the CLR do not contradict with the locality and 
the translational invariance on lattice. Yet in the naive continuum limit, CLR cannot be 
distinguished with LR. Because of these features we consider the CLR may be an important 
key ingredient for further constructing super symmetric lattice models beyond the quantum 
mechanics model discussed in this paper. 

Using a difference operator and field products which satisfy the CLR, we can easily 
construct exact half-supersymmetric action with superfield formalism. We can also find 
the Nicolai map without "surface term". Since each term of our superspace action is 
independently invariant under supersymmetry, we can introduce deformation parameters 
with no difficulty. This enable us to utilize a localization technique to compute, for instance, 
the Witten index of our supersymmetric quantum mechanics models, which reproduces the 
known results in the continuum. 

In order to construct supersymmetric field theoretical model, we have to adapt spinor 
index for higher dimensional CLR. In this paper, as a first step, we consider two-dimensional 
CLR without spinor index, in which we find a difficulty that the difference operator has a 
particular form unsuitable for the kinetic term. We hope to report a further analysis which 
properly includes spinor index elsewhere. 
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A. General solution of the CLR and the properties of smeared symmetric 
product 

A general difference operator is written as a function of the difference of site index 

A mn = A(m-n), (A.l) 

by the translational invariance. The momentum representation (if-expression) is written 
ELS cl holomorphic function, 

oo 

A(w) = w n A{n), (A.2) 

n=— oo 

where w = exp(ipa) in terms of the momentum p dual to the coordinate na. This is 
holomorphic in the domain D which is an annulus including a unit circle around the origin. 
From (2.3), the constraint 

A(l) = (A.3) 

is imposed on the difference operator. For the product coefficient between lattice fields, we 
can also write in the translational invariant form 

M lmn = M(£-m,e-n). (A A) 

The momentum representation is written as a two-variables holomorphic function, 



DC 



M(w,z)= w m z n M(m,n) = M{z,w), (A.5) 



m,n=—oo 



where w = exp(ipa), z = exp(iqa) and the holomorphic domain is T>®T>. Then, our cyclic 
Leibniz rule ( p.21| ) can be expressed as 

M(w, z)A(-) + M(z, —)A(w) + M(—,w)A(z) = 0. (A.6) 

wz wz wz 

The general solution for this equation is 

MO, z)A(-) = f(w, z) + f(z, 1) + f(wz, — ) + f(z, w) + f(w, 1) + f(wz, — ) 
wz wz wz 

- f(z, — ) - /(— , 1) - f(-,w) - f(w, — ) - /(— , 1) - /(-, z), 

wz wz w wz wz z 

where f(w,z) is an arbitrary two-variable holomorphic function. By taking f(w,z) = 
F(w) + F(z), the special solution is found 

M(w, z)A( — ) = F(w) - F(-) + F(z) - F(-) - 2F(—) + 2F(wz). (A.7) 

wz w z wz 

2 

If we choose F(w) = j^, then 

MKz)A(i) = 2^ + ^ + ^- 1 + 2W 1 u;z- W 1 

v ' 1 y wz' 6 2 K ' 
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This corresponds to the momentum representation of ( 2.25 ). 

Note that we cannot impose an associative law on our symmetric product. Namely, 
the product QA.4 ) cannot satisfy the associative law 



W^,x)} = {{Wk}. (A.9) 



If we imposed ( A.9 ), then 

M £mn M jek = ^ M jm£ M £nk (A.10) 
i t 

or in tu-expression, 

M(w, z)M(wz, x) = M(z, x)M(xz, w) (All) 
would hold which implied M(w,z) = 1 after local field redefinition 0], while the constant 



M was excluded from our general solution (A/7) . 



Another property of the symmetric product ( |A.4| ) may come out for a case where there 
exists a unit element I which satisfies 

{</>, /} = {/, <p} = 4>. (A.12) 

Since I n can be taken as a unit constant by a local field redefinition, combining with the 
translational invariance, the above implies 

M(w, 1) = M(l, w) = l (A.13) 



for any w. Using the solution (A. 7), we can find an example satisfying (A.13) 



M(w, z) = 1 + —^—(w -l)(z- l)(w + z-2- 2wz). (A. 14) 

12wz 

B. Multi-body products and the CLR 

In this appendix we describe how multi-body smeared symmetric products satisfying CLR 
can be constructed from 2-body ones. 

B.l Construction 

A^-body smeared symmetric product {i^ 1 -*, (j)^ 2 \ • • • , (j)^} is defined as 

ni, ■■■,n jv 

with the coefficient M n;ni ... njv whose last indices are totally symmetric. The CLR relation 
for this can be expressed as 

(A<M</v,</>}) = (B.2) 

for a single field eft, or 

(A0 c «,{^ 2 ),...,^ +1 )}) =0 (B.3) 

{c(i) : cyclic perm.} 
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for N + 1 fields (j) 1 (i = 1,2, • • • , iV + 1), where the sum is taken over all possible cyclic 
permutations {c(i) \ i = 1, • • • , N + 1}. 

iV-body symmetric product can be constructed from 2-body product by the following 
steps: 

1. Write down a tree diagram (without loops) with N + 1 external lines consists only of 
trivalent vertices. (See Fig. |l| for aniV = 4 example.) 




Figure 1: Example of a tree diagram from which one can construct 4-body product 



2. Choose one of the external lines and assign it "0". Then number the rest of external 
lines from 1 to N in arbitrary order. (See Fig. ^ for counter-clockwise order examples) 

3. For each vertex mark a dot on a line, among three lines attached to the vertex, 
which is connected to the portion containing "0"-line. (See Fig. ^) Now this dia- 
gram expresses one of 4-body products composed of three 2-body products. Each 
vertex corresponds to 2-body product coefficient M nm i and a dot distinguishes first 
index. Indices of internal lines (lines connecting vertices) are summed. Thus, for 
example, Fig. |2|(a) corresponds to Y2k i ^nn^k^kn^i^inzn^ which gives a product of 
fields {0( 2 ), (4) }}}. 




Figure 2: Constructing 4-body product made of 2-body products 



4. For all possible choices of "0"-line and any permutatios of "1" through "iV", repeat 
the steps 2 and 3. Then sum up all combined products corresponding to these 
diagrams with normalization factor 1/(N + 1)!. For N = 4 case, all possible choices 
of "0"-line obtained from Fig. p] are shown in Fig. ^, and the final expression is 

M n]nin2n3n4 = — (M nnva) kMk nv(2) iMi nv(3)nriU) + Mk npW iMi np ^ np ^M n kn pW 

' perm. k,l 
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+ -^nfc/-^fcnp( 1 )np(2)^np(3)n p (4) + ^nn p{1) k^ln p (^n p{:i) ^kln pW 
+ M ln p(1) n p(2) M k i nm M nknpi4) j . (B.4) 
This gives 4-body field product as 

{^\^,^\^} = i y: ({^ (i)) ,{^ (2)) ,{^ (3)) ,^ (4)) }}} 

' perm. 

, {</X 2 » , </> (3)) }} , </> (p(4)) } + , (p(2)) } , {</> (p(3)) , </> (4)) }} 

+{0W 1 )),{{0W 2 )),0^ 3 ))},0W 4 ))}} + {{{^W 1 )),^ 2 ))},^ 3 ))},^ 4 ))}) .(B.5) 



Thus we obtain one of possible definitions of multi-body symmetric product for each tree 
diagram we start from. Although tree diagram for N = 3 or 4 case is unique, there are 
a finite number of independent tree diagrams for higher N. Each of these diagrams gives 
independent product on a lattice, while they may have the same continuum limit. 

Note that a class of multi-body symmetric products can be alternatively obtained by 
a simple equation 

1 - 

{0, <t>, ■ ■ ■ A} = jy-q— j" (<t>), Pk(<f>)}- (B.6) 

k=0 

Here the field P n {4>) is defined in a recursive way using 2-body product 

Pn{4>) = {P„_i(0),0}, Pi{4>) = ^ Po(<P) = I (B.7) 
where I is an identity field with respect to the 2-body product (eft = {/,</>}). 
B.2 Proof of multi-body CLR 

Here we show that the above constructed products satisfy the multi-body CLR. A graphical 
representation is convenient for this purpose. As in the previous subsection, coefficient of 
2-body product M nrn x is expressed by the trivalent vertex with a dot on the first index. 
Difference operator A nm is now expressed by a line segment with an arrow: 

/ 

M nm i n A nm n < m (B.8) 

Anti-symmetry of A (= — A T ) is represented as a graphical relation: 

" < rn = - n > m ^ g . 

Then CLR for 2-body product X^cyclic (^^ c > {4> c<y2 \ (f> c ^}) = is graphically repre- 
sented by Fig. |||. 
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Figure 3: CLR for 2-body products 



Let us begin by looking at N = 4 example. What we want to show is the relation 



: cyclic perms. 



(B.10) 



Substituting the definition of 4-body product (EL5) into the left-hand side of (|B.1C| ), we 
have 600 terms in total (5! kinds of terms times 5 cyclic permutations). These can be 
rearranged into 5! groups, one of which made of five graphs shown in Fig. || and the others 
are its permutations with respect to external lines. Therefor it is enough to argue the 
cancellations of graphs in one group. Actually contributions of five graphs in Fig. [| add up 
to zero by the CLR for 2-body product (Fig. |3]) and the anti-symmetry of A. For example, 




Figure 4: Vanishing set of graphs in CLR for 4-body products 




first two graphs in Fig. [| add up to become the graph Fig. ||(a) with minus sign by the 2- 
body CLR, while last two graphs in Fig. || do become the graph Fig. ||(b) again with minus 



sign. Remaining Fig. ||(c), however, can be rewritten into a sum of the graphs Fig. |5|(c) 
and (d) with minus signs. Therefor they cancel out altogether due to anti-symmetry of A. 

We can generalize the above arguments to arbitrary N case by induction. Suppose we 
already have the M-body CLR for M less than N. Consider a tree graph with N+l = K+L 
external lines (See left hand side of Fig. |^). We can find an internal line which separates 
the graph into two parts each of which has K and L external lines respectively (See right 
hand side of Fig. ||), where 1 < K < N and 1 < L < N. To examine the CLR, we attach 




Figure 6: Factorization into products with smaller numbers 



A at one of external lines and sum up each contributions. Then the contributions of the 
graphs attached A at an external line with number 1 though K are turned out to be that 
of Fig. ^(a) with minus sign by the K-ho&y CLR, while those with number K + 1 through 
K + L be that of Fig. 0(b) with minus sign by the -L-body CLR. As was the case of N = 4, 
Fig. ^(a) and (b) cancel each other by the anti-symmetry of A. Thus iV-body CLR is 
satisfied. 




Figure 7: Canceling graphs: (a)+(b)=0 due to A + A T = 
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